Computational methods meet in vitro techniques: A case study on fusaric acid and its possible detoxification through cytochrome P450 enzymes.

Mycotoxins are known environmental pollutants that may contaminate food and feed chains. Some mycotoxins are regulated in many countries to limit the trading of contaminated and harmful commodities. However, the so-called emerging mycotoxins are poorly understood and need to be investigated further. Fusaric acid is an emerging mycotoxin, noxious to plants and animals, but is known to be less toxic to plants when hydroxylated. The detoxification routes effective in animals have not been elucidated yet. In this context, this study integrated in silico and in vitro techniques to discover potential bioremediation routes to turn fusaric acid to its less toxic metabolites. The toxicodynamics of these forms in humans have also been addressed. An in silico screening process, followed by molecular docking and dynamics studies, identified CYP199A4 from the bacterium Rhodopseudomonas palustris HaA2 as a potential fusaric acid biotransforming enzyme. Its activity was confirmed in vitro. However, the effect of hydroxylation seemed to have a limited impact on the modelled toxicodynamics against human targets. This study represents a starting point to develop a hybrid in silico/in vitro pipeline to find bioremediation agents for other food, feed and environmental contaminants.

Free fatty acid receptors beyond fatty acids: A computational journey to explore peptides as possible binders of GPR120.

Free fatty acids receptors, with members among G protein-coupled receptors (GPCRs), are crucial for biological signaling, including the perception of the so called "fatty taste". In recent years, GPR120, a protein belonging to the GPCR family, drew attention as an interesting pharmacological target to cope with obesity, satiety and diabetes. Apart from long chain fatty acids, which are GPR120 natural agonists, other synthetic molecules were identified as agonists expanding the chemical space of GPR120's ligands. In this scenario, we unveiled peptides as possible GPR120 binders toward a better understanding of this multifaceted and relevant target. This study analyzed a virtual library collecting 531 441 low-polar hexapeptides, providing mechanistic insights on the GPR120 activation and further extending the possible chemical space of GPR120 agonists. The computational pipeline started with a narrow filtering of hexapeptides based on their chemical similarity with known GPR120 agonists. The best hits were tested through docking studies, molecular dynamics and umbrella sampling simulations, which pointed to G[I,L]FGGG as a promising GPR120 agonist sequence. The presence of both peptides in food-related proteins was thoroughly assessed, revealing they may occur in mushrooms, food-grade bacteria and rice. Simulations on the counterparts with D-amino acids were also performed. Umbrella sampling simulations described that GdIFGGG may have a better interaction compared to its all-L counterpart (-13 kCal/mol ΔG and -6 kCal/mol ΔG, respectively). Overall, we obtained a predictive model to better understand the underpinning mechanism of GPR120-hexapeptides interaction, hierarchizing novel potential agonist peptides for further analysis and describing promising food sources worth of further dedicated investigations.

Virtual display of targets: A new level to rise the current understanding of ochratoxin A toxicity from a molecular standpoint.

Ochratoxin A (OTA) is a mycotoxin spread worldwide contaminating several food and feed commodities and rising concerns for humans and animals. OTA toxicity has been thoroughly assessed over the last 60 years revealing a variety of adverse effects, including nephrotoxicity, hepatotoxicity and possible carcinogenicity. However, the underpinning mechanisms of action have yet to be completely displayed and understood. In this framework, we applied a virtual pipeline based on molecular docking, dynamics and umbrella simulations to display new OTA potential targets. The results collected consistently identified OGFOD1, a key player in protein translation, as possibly inhibited by OTA and its 2'R diastereomer. This is consistent with the current knowledge of OTA's molecular toxicology and may fill some gaps from a mechanistic standpoint. This could pave the way for further dedicated analysis focusing their attention on the OTA-OGFOD1 interaction, expanding the current understanding of OTA toxicity at a molecular level.

Occurrence of mycotoxins in meat alternatives: Dietary exposure, potential health risks, and burden of disease.

This study aimed to present the occurrence of sixteen mycotoxins in 105 meat alternatives based on wheat, legumes, and vegetables from Italy. The targeted mycotoxins were aflatoxins (AFB1, AFB2, AFG1, AFG2), fumonisins B1 and B2 (FB1, FB2), alternariol (AOH), alternariol monomethyl ether (AME), tentoxin (TEN), ochratoxin A (OTA), zearalenone (ZEN), T-2/HT-2 toxin, deoxynivalenol (DON), enniatin B (ENNB), and beauvericin (BEA). The occurrence of mycotoxins was between 0% (AFB2) - 97.4% (ENNB). Mycotoxin co-occurrence varied from binary combinations up to mixtures of twelve. To assess the dietary exposure and potential health risks we simulated the replacement of meat consumption for Italian consumers with meat alternatives. The cumulative exposure to Alternaria mycotoxins and trichothecenes indicated a potential health risk while the exposure to aflatoxins and ochratoxin A indicated a potential health concern related to liver and renal cancer in the model scenario. Moreover, we estimated the risk of liver cancer from exposure to AFB1 and quantified the potential burden using Disability-Adjusted Life Years (DALYs). Luckily, the potential risk of liver cancer was low between 0 and 0.05/100,000 individuals with an associated burden of disease of 0.83 DALYs/100,000 individuals. Taking into consideration the presence of meat alternatives on the food market and the ongoing shift towards plant-based diets there is a need for continuous monitoring to keep the occurrence at safe levels. More attention is needed from the regulatory side for policymakers to consider the legislations of mycotoxins in meat alternatives.

Investigating Metabolic Plant Response toward Deoxynivalenol Accumulation in Four Winter Cereals.

Deoxynivalenol (DON) is a phytotoxic agent supporting the spread of fungal diseases in cereals worldwide, i.e., fusarium head blight. It is known that DON accumulation may elicit changes in plant secondary metabolites in response to pathogen attack. This study maps the changes in selected secondary metabolite classes upon DON contamination occurring in fifteen Triticum spp. genotypes, among them emmer, spelt, and soft wheat, and 2 tritordeum varieties, cultivated in two different sites and over two harvest years. The main phenolic classes (i.e., alkylresorcinols, soluble, and cell-wall bound phenolic acids) were targeted analyzed, while changes in the lipidome signature were collected through untargeted HRMS experiments. The results, obtained across multiple Triticum species and in open fields, confirmed the modulation of first-line biological pathways already described in previous studies involving single cereal species or a limited germplasm, thus reinforcing the involvement of nonspecific chemical defenses in the plant response to pathogen attack.

Mycotoxin occurrence in kernels and straws of wheat, barley, and tritordeum.

Thirty-two varieties of common and durum wheat, hordeum, barley, and tritordeum collected over two harvesting years (2020 and 2021) were investigated for the presence of multiple Fusarium-related mycotoxins in asymptomatic plants. DON, 3-AcDON, 15-AcDON, T-2, HT-2, and ZEN together with the emerging mycotoxin ENN B and the major modified form of DON, namely DON3Glc, were quantified by means of UHPLC-MS/MS. Overall, DON and ENN B were the most frequently detected mycotoxins, albeit large inter-year variability was observed and related to different climate and weather conditions. Straws had higher mycotoxin contents than kernels and regarding DON occurrence tritordeum was found to be the most contaminated group on average for both harvesting years, while barley was the less contaminated one. Emerging mycotoxin ENN B showed comparable contents in kernels compared to straw, with a ratio close to 1 for tritordeum and barley. Regarding the occurrence of the other evaluated mycotoxins, T-2 and HT-2 toxins have been spotted in a few tritordeum samples, while ZEN has been frequently found only in straw from the harvesting year 2020. The data collected confirms the occurrence of multiple Fusarium mycotoxins in straws also from asymptomatic plants, highlighting concerns related to feed safety and animal health. The susceptibility of Tritordeum, hereby reported for the first time, suggests that careful measures in terms of monitoring, breeding, and cultural choices should be applied when dealing with his emerging crop.

A mechanistic toxicology study to grasp the mechanics of zearalenone estrogenicity: Spotlighting aromatase and the effects of its genetic variability.

Zearalenone (ZEN) is a mycoestrogen produced by Fusarium fungi contaminating cereals and in grain-based products threatening human and animal health due to its endocrine disrupting effects. Germane to the mechanisms of action, ZEN may activate the estrogen receptors and inhibit the estrogens-producing enzyme aromatase (CYP19A1). Both show single nucleotide variants (SNVs) among humans associated with a diverse susceptibility of being activated or inhibited. These variations might modify the endocrine disrupting action of ZEN, requiring dedicated studies to improve its toxicological understanding. This work focused on human aromatase investigating via 3D molecular modelling whether some of the SNVs reported so far (n = 434) may affect the inhibitory potential of ZEN. It has been also calculated the inhibition capability of α-zearalenol, the most prominent and estrogenically potent phase I metabolite of ZEN, toward those aromatase variants with an expected diverse sensitivity of being inhibited by ZEN. The study: i) described SNVs likely associated with a different susceptibility to ZEN and α-zearalenol inhibition - like T310S that is likely more susceptible to inhibition, or D309G and S478F that are possibly inactive variants; ii) proofed the possible existence of inter-individual susceptibility to ZEN; iii) prioritized aromatase variants for future investigations toward a better comprehension of ZEN xenoestrogenicity at an individual level.

Untargeted Metabolomics of Meat Digests: Its Potential to Differentiate Pork Depending on the Feeding Regimen.

Meat quality seems to be influenced by the dietary regimes applied for animal feeding. Several research studies are aimed at improving meat quality, preserving it from oxidative processes, by the incorporation of antioxidant components in animal feeding. The main part of these studies evaluates meat quality, determining different parameters directly on meat, while few research studies take into account what may happen after meat ingestion. To address this topic, in this study, an in vitro gastrointestinal digestion protocol was applied to two different pork muscles, longissimus dorsi and rectus femoris, obtained from pigs fed with different diets. In detail, two groups of 12 animals each were subjected to either a conventional diet or a supplemented diet with extruded linseeds as a source of omega-3 fatty acids and plant extracts as a source of phenolics antioxidant compounds. The digested meat was subjected to an untargeted metabolomics approach. Several metabolites deriving from lipid and protein digestion were detected. Our untargeted approach allowed for discriminating the two different meat cuts, based on their metabolomic profiles. Nonetheless, multivariate statistics allowed clearly discriminating between samples obtained from different animal diets. In particular, the inclusion of linseeds and polyphenols in the animal diet led to a decrease in metabolites generated from oxidative degradation reactions, in comparison to the conventional diet group. In the latter, fatty acyls, fatty aldehydes and oxylipins, as well as cholesterol and vitamin D3 precursors and derivatives, could be highlighted.

The bitter side of toxicity: A big data analysis spotted the interaction between trichothecenes and bitter receptors.

The bitter taste perception evolved in human and animals to rapidly perceive and avoid potential toxic compounds. This is mediated by taste receptors type 2 (TAS2R), expressed in various tissues, which recently proved to be involved in roles beyond the bitter perception itself. With this study, the interaction between food-related toxic compounds and TAS2R46 has been investigated via computational approaches, starting with a virtual screening and moving to molecular docking and dynamics simulations. The virtual screening analysis identified trichothecolone and the trichothecenes class it belongs to, which includes mycotoxins widespread in several commodities raising food safety concerns, as possible TAS2R46 binders. Molecular docking and dynamics simulations were performed to further explore the trichotecenes-TAS2R46 interaction. The results indicated that deoxynivalenol and its 15-acetylated derivative could activate TAS2R46. Eventually, this study provided initial evidence supporting the involvement of TAS2R46 in the underpinning mechanisms of deoxynivalenol action highlighting the need of digging into the involvement of TAS2R46 and TAS2Rs in the adverse effects of deoxynivalenol and congeners.

Occurrence of Ergot Alkaloids in Major and Minor Cereals from Northern Italy: A Three Harvesting Years Scenario.

Ergot alkaloids (EAs), mycotoxins produced mainly by fungi of the Claviceps genus, have been frequently reported in rye, while their increasingly frequent occurrence in other cereals is likely related to weather conditions, with the incidence of ergot sclerotia in winter grains being related to heavy rainfall and moist soils at critical periods. However, compared to other regulated mycotoxins, data about the prevalence and occurrence of EAs in major and minor cereals harvested in the Mediterranean growing areas are still scant. In this regard, the current study reported the occurrence of EAs in 18 genotypes of winter cereals harvested over 3 years from an experimental field located in North Italy which were analyzed by HPLC-MS/MS. Results indicate a widespread occurrence of all the major EAs in all the considered cereal crops, especially under supportive meteorological conditions. EA contamination was dependent on the harvest year (p < 0.0001) which was particularly high in 2020 for all the considered species. The results also demonstrated a large co-occurrence of EAs with 98 cereal samples out of 162 contaminated with at least one of the 12 EAs (60% positive samples) in the range LOD: 15,389 µg/kg (median value: 2.32 µg/kg), expressed as the sum of the EAs. Rye was confirmed to be the crop more susceptible to the fungal infection (EAs content up to 4,302 µg/kg). To the best of our knowledge, we have reported the accumulation of EAs in tritordeum (LOD: 15,389 µg/kg) and in emmer (LOD: 1.9 µg/kg) for the first time.

The Occurrence of Non-Regulated Mycotoxins in Foods: A Systematic Review.

The aim of this systematic review is to provide an update on the occurrence and co-occurrence of selected non-regulated mycotoxins and provide an overview of current regulations. Fifteen non-regulated mycotoxins were found in 19 food categories worldwide. On top of that, 38 different combinations of non-regulated mycotoxins were found, with mixtures varying from binary combinations up to 12 mycotoxins. Taking into consideration the amount of evidence regarding the prevalence and co-occurrence of non-regulated mycotoxins, future steps should be taken considering continuous monitoring, scientific exchange, and generation of high-quality data. To enhance data quality, guidelines outlining the minimum quality criteria for both occurrence data and metadata are needed. By doing so, we can effectively address concerns related to the toxicity of non-regulated mycotoxins. Furthermore, obtaining more data concerning the co-occurrence of both regulated and non-regulated mycotoxins could aid in supporting multiple chemical risk assessment methodologies. Implementing these steps could bolster food safety measures, promote evidence-based regulations, and ultimately safeguard public health from the potential adverse effects of non-regulated mycotoxins.

An in silico insight on the mechanistic aspects of gelsenicine toxicity: A reverse screening study pointing to the possible involvement of acetylcholine binding receptor.

Gelsedine-type alkaloids are highly toxic plant secondary metabolites produced by shrubs belonging to the Gelsemium genus. Gelsenicine is one of the most concerning gelsedine-type alkaloids with a lethal dose lower than 1 mg/Kg in mice. Several reported episodes of poisoning in livestock and fatality cases in humans due to the usage of Gelsemium plants extracts were reported. Also, gelsedine-type alkaloids were found in honey constituting a potential food safety issue. However, their toxicological understanding is scarce and the molecular mechanism underpinning their toxicity needs further investigations. In this context, an in silico approach based on reverse screening, docking and molecular dynamics successfully identified a possible gelsenicine biological target shedding light on its toxicodynamics. In line with the available crystallographic data, it emerged gelsenicine could target the acetylcholine binding protein possibly acting as a partial agonist against α7 nicotinic acetylcholine receptor (AChR). Overall, these results agreed with evidence previously reported and prioritized AChR for further dedicated analysis.

Untargeted metabolomics liquid chromatography-high resolution mass spectrometry approach for the geographical origin assessment of Italian dehydrated apples.

Geographical provenience is nowadays a relevant aspect of the authenticity and the quality of many food commodities. Dehydrated apple cubes/slices represent an ingredient commonly used by food companies for bakery products. However, this apple-based matrix is not so known and studied from an analytical point of view. In the present work, seven compounds were identified as key molecules to distinguish between Italian and non-Italian samples, through an untargeted ultrahigh-pressure liquid chromatography-high-resolution mass spectrometry (UHPLC-HRMS) approach. This methodology was merged with multivariate statistical analysis, and the principal features were studied and identified considering several identification steps. Samples from 2020 and 2021 harvesting campaigns, with partial and total dehydration rates, with or without peel, and from different apple varieties were considered for the study, for a total of 91 samples. Afterward, the same analysis protocol was applied to an external set (n = 12 samples), included in the statistical models, searching for the key compounds identified in the training set. Interesting and significant results underlined the potentiality of the UHPLC-HRMS technology as a confirmatory strategy for the geographical origin assessment of dehydrated apple commodities.

QuEChERS method combined to liquid chromatography high-resolution mass spectrometry for the accurate and sensitive simultaneous determination of pyrrolizidine and tropane alkaloids in cereals and spices.

Within the last decades, in the EU, there has been an increasing interest in toxic plant alkaloids as food contaminants, especially after the continuous and growing consumption of plant-based foods compared with food of animal origin. In this regard, the once neglected presence of these tropane alkaloids (TAs) and pyrrolizidine alkaloids (PAs) has recently been reconsidered by the European Food Safety Authority, highlighting the lack of data and the need to develop risk assessment strategies. For this reason, the emphasis has been placed on detecting their occurrence in food through the development of accurate and sensitive analytical methods to achieve the determination of these compounds. The present study aims to elaborate and validate an analytical method based on QuEChERS sample preparation approach, exploiting the UHPLC coupled to the HRMS to simultaneously identify and quantify 21 PAs and two TAs in cereals and spices. For TAs, the obtained limit of detection (LOD) is 0.1 µg·kg-1 and the limit of quantification (LOQ) is 0.4 µg·kg-1 , while for PAs, the LODs values ranging between 0.2 to 0.3 µg·kg-1 and the LOQ, between 0.4 and 0.8 µg·kg-1 , ensuring compliance with the recently established European Regulations. Several commercial samples were analysed to further verify the applicability of this comprehensive analytical approach.

Food processing contaminants: Dietary exposure to 3-MCPD and glycidol and associated burden of disease for Italian consumers.

In this study we assessed Italian consumers' dietary exposure to 3-MCPD and glycidol followed by risk characterization, potential cancer risk and the associated burden of disease. Consumption data was retrieved from the most recent Italian Food Consumption Survey (2017-2020), while contamination data was obtained from the European Food Safety Authority. The level of risk due to exposure to 3-MCPD was negligible, below the tolerable daily intake (TDI), except for high consumption of infant formulas. For infants, the intake level was higher than the TDI (139-141% of TDI), indicating a potential health risk. Exposure to glycidol indicated a health concern for infants, toddlers, other children, and adolescents consuming infant formulas, plain cakes, chocolate spreads, processed cereals, biscuits, rusks, and cookies (margin of exposure (MOE) < 25,000). The risk of cancer due to exposure to glycidol was estimated and the overall health impact was quantified in Disability-Adjusted Life Years (DALYs). The risk of cancer due to chronic dietary exposure to glycidol was estimated at 0.08-0.52 cancer cases/year/100,000 individuals depending on the life stage and dietary habits in Italy. The burden of disease quantified in DALYs varied from 0.7 to 5.37 DALYs/year/100,000 individuals. It is crucial to continuously gather consumption and occurrence data for glycidol over time to track patterns, assess potential health risks, identify exposure sources, and develop effective countermeasures, as long-term exposure to chemical contaminants can lead to an increased risk for human health. This data is critical for protecting public health and reducing the likelihood of cancer and other health issues related to glycidol exposure.

A geographical origin assessment of Italian hazelnuts: Gas chromatography-ion mobility spectrometry coupled with multivariate statistical analysis and data fusion approach.

Hazelnut is a commodity that has gained interest in the food science community concerning its authenticity. The quality of the Italian hazelnuts is guaranteed by Protected Designation of Origin and Protected Geographical Indication certificates. However, due to their modest availability and the high price, fraudulent producers/suppliers blend, or even substitute, Italian hazelnuts with others from different countries, having a lower price, and often a lower quality. To contrast or prevent these illegal activities, the present work investigated the application of the Gas Chromatography-Ion mobility spectrometry (GC-IMS) technique on the hazelnut chain (fresh, roasted, and paste of hazelnuts). The raw data obtained were handled and elaborated using two different ways, software for statistical analysis, and a programming language. In both cases, Principal Component Analysis and Partial Least Squares-Discriminant Analysis models were exploited, to study how the Volatile Organic Profiles of Italian, Turkish, Georgian, and Azerbaijani products differ. A prediction set was extrapolated from the training set, for a preliminary models' evaluation, then an external validation set, containing blended samples, was analysed. Both approaches highlighted an interesting class separation and good model parameters (accuracy, precision, sensitivity, specificity, F1-score). Moreover, a data fusion approach with a complementary methodology, sensory analysis, was achieved, to estimate the performance enhancement of the statistical models, considering more discriminant variables and integrating at the same time further information correlated to quality aspects. GC-IMS could be a key player as a rapid, direct, cost-effective strategy to face authenticity issues regarding the hazelnut chain.

Implementing the Use of Collision Cross Section Database for Phycotoxin Screening Analysis.

The increased consumption of blue-green algae (BGA)-based dietary supplements has raised concern about their food safety, especially about cyanotoxin presence. The hyphenation of liquid chromatography with ion mobility mass spectrometry represents a relevant tool to screen several compounds in a large variety of food matrices. In this work, ultrahigh-performance liquid chromatography coupled to traveling wave ion mobility spectrometry/quadrupole time-of-flight mass spectrometry (UHPLC-TWIMS-QTOF) was employed to establish the first comprehensive TWIMS-derived collision cross section database (TWCCSN2) for phycotoxins. The database included 20 cyanotoxins and 1 marine toxin. Accurate m/z, retention times, and TWCCSN2 values were obtained for 81 adducts in positive and negative electrospray (ESI+/ESI-) modes. Reproducibility and robustness of the TWCCSN2 measurements were determined to be independent of the matrix. A screening was carried out on 19 commercial BGA dietary supplements of different composition. Cyanotoxins were confidently identified in five samples based on retention time, m/z, and TWCCSN2.

Multi-mycotoxin determination in plant-based meat alternatives and exposure assessment.

The aim of this study was to fill in the gap regarding the occurrence of mycotoxins in plant-based meat alternatives. Hence, a multi-mycotoxin method (aflatoxins, ochratoxin A, fumonisins, zearalenone, and mycotoxins from the Alternaria alternata genera) was developed followed by an exposure assessment for the Italian consumers' exposure to mycotoxins. A total of 13 meat alternatives samples based on soy, pea, chickpea, lupin, and seitan were analysed. With the exception of seitan, all of the remaining samples were contaminated with one mycotoxin or mixtures of up to seven mycotoxins. The level of contamination was as low as 0.2 µg/kg alternariol methyl ether and as high as 66.9 µg/kg fumonisin B1. To analyse the exposure to mycotoxins due to plant-based meat alternatives consumption we used the consumption meat data from the Food and Agriculture Organization for Italian adult consumers and simulated a full replacement of meat with plant-based meat alternatives. Based on our model, consumption of plant-based meat alternatives led to a non-tolerable exposure to alternariol (hazard index (HI) > 1) in pea-based burger and soy + wheat-based steak, while samples contaminated with aflatoxins, respectively ochratoxin A, indicated a health concern related to liver and renal cancer (margin of exposure (MOE) < 10,000). This is the first study that presents the co-occurrence of mycotoxins in multiple plant-based meat alternatives. Moreover, these results indicate that there is a need for policymakers to consider the regulation of mycotoxins in plant-based meat alternatives in order to ensure consumers' safety.

A computational study on the biotransformation of alkenylbenzenes by a selection of CYPs: Reflections on their possible bioactivation.

Alkenylbenzenes are aromatic compounds found in several vegetable foods that can cause genotoxicity upon bioactivation by members of the cytochrome P450 (CYP) family, forming 1'-hydroxy metabolites. These intermediates act as proximate carcinogens and can be further converted into reactive 1'-sulfooxy metabolites, which are the ultimate carcinogens responsible for genotoxicity. Safrole, a member of this class, has been banned as a food or feed additive in many countries based on its genotoxicity and carcinogenicity. However, it can still enter the food and feed chain. There is limited information about the toxicity of other alkenylbenzenes that may be present in safrole-containing foods, such as myristicin, apiole, and dillapiole. In vitro studies showed safrole as mainly bioactivated by CYP2A6 to form its proximate carcinogen, while for myristicin this is mainly done by CYP1A1. However, it is not known whether CYP1A1 and CYP2A6 can activate apiole and dillapiole. The present study uses an in silico pipeline to investigate this knowledge gap and determine whether CYP1A1 and CYP2A6 may play a role in the bioactivation of these alkenylbenzenes. The study found that the bioactivation of apiole and dillapiole by CYP1A1 and CYP2A6 is limited, possibly indicating that these compounds may have limited toxicity, while describing a possible role of CYP1A1 in the bioactivation of safrole. The study expands the current understanding of safrole toxicity and bioactivation and helps understand the mechanisms of CYPs involved in the bioactivation of alkenylbenzenes. This information is essential for a more informed analysis of alkenylbenzenes toxicity and risk assessment.

A Computational Inter-Species Study on Safrole Phase I Metabolism-Dependent Bioactivation: A Mechanistic Insight into the Study of Possible Differences among Species.

Safrole, a 162.2 Da natural compound belonging to the alkenylbenzenes class, is classified as a possible carcinogen to humans by IARC (group IIB) and has proven to be genotoxic and carcinogenic to rodents. Despite its use as a food or feed additive, it is forbidden in many countries due to its documented toxicity; yet, it is still broadly present within food and feed and is particularly abundant in spices, herbs and essential oils. Specifically, safrole may exert its toxicity upon bioactivation to its proximate carcinogen 1'-hydroxy-safrole via specific members of the cytochrome P450 protein family with a certain inter/intra-species variability. To investigate this variability, an in-silico workflow based on molecular modelling, docking and molecular dynamics has been successfully applied. This work highlighted the mechanistic basis underpinning differences among humans, cats, chickens, goats, sheep, dogs, mice, pigs, rats and rabbits. The chosen metric to estimate the likeliness of formation of 1'-hydroxy-safrole by the species-specific cytochrome P450 under investigation allowed for the provision of a knowledge-based ground to rationally design and prioritise further experiments and deepen the current understanding of alkenylbenzenes bioactivation and CYPs mechanics. Both are crucial for a more informed framework of analysis for safrole toxicity.